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Coupling  Problems  with  High  Birefringence  Fibers 
Chin>Lin  Chen 

School  of  Electrical  Engineering 
Purdue  University 
West  Lafayette,  IN  47007 

Abstract 

A  program  to  investigate  the  coupling  problems  of  high  birefringence  fibers  was 
initiated  «,in  1983  for  Naval  Research  Laboratory  under  Contract  N00014-83-K-2029. 
The  technical  monitor  was  Dr.  W.  K.  Burns  of  Naval  Research  Laboratory.  The  work 
performed  under  this  contract  is  summarized  in  this  report. 

A  fiber  is  birefringent  if  it  has  a  non-circular  core  and/or  cladding,  or  is  subjected 
to  anisotropic  stress.  Currently,  the  loss  of  fibers  with  embedded  anisotropic  stress  is 
much  lower  than  that  with  non-circular  core  or  cladding  and  therefore  these  fibers  are 
of  interest.  Central  to  any  coupling  problem  is  the  distribution  of  fields  in  the  fibers. 
To  date,  very  little  is  known  relative  to  the  fields  of  polarization  modes  in  biaxial 
fibers.  This  is  the  main  thrust  of  our  study.  We  have  used  a  perturbation  method  to 
study  the  fields  and  dispersion  of  polarization  modes  in  biaxial  fibers.  In  our  approach, 
the  inclusion  or  exclusion  of  a  term  is  determined  by  the  differential  equation,  the 
boundary  conditions  and  the  symmetry  of  the  problem,  and  does  not  rely  on  the  prior 
knowledge  of  the  fields  in  the  fibers.  The  dependence  of  the  fiber  birefringence  on  the 
index  differential  and  index  anisotropy  has  been  studied.  Now  that  the  fields  in  high 
birefringent  fiber  are  known,  the  excitation  of  high  birefringent  fibers  can  be  analyzed. 
Methods  for  studying  these  excitation  problems  are  indicated  and  some  of  the  results 
are  presented. 
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1.  Introduction 


Since  low-loss  fibers  were  realized  some  two  decades  ago,  interest  on  fibers  has 
been  evolved  from  multimode  fibers  to  single  mode  fibers.  In  reality,  all  conventional 
single  mode  fibers  support  two  orthogonal  polarization  modes.  As  these  polarization 
modes  arc  nearly  degenerate,  fields  in  these  fibers  are  easily  converted  from  one  polari¬ 
zation  mode  to  the  another  if  there  is  slightly  disturbance  or  perturbation.  For  sophis¬ 
ticated  systems  like  coherent  communication  systems,  sensory  systems  based  on  inter¬ 
ferometric  principles,  or  systems  utilizing  polarization-dependent  components,  these 
changes  in  the  state  of  polarization  would  lead  to  signal  fading  and  noise.  The  polari¬ 
zation  mode  conversion,  being  the  result  of  mode  coupling,  can  be  greatly  reduced  if 
the  polarization  mode  degeneracy  is  removed  [1,2].  With  polarization  mode  degeneracy 
removed,  fibers  become  highly  birefringent  since  their  propagation  velocities  are  quite 
differenl.  If  the  birefringence  is  sufTiciently  large,  and  the  fiber  core  is  sufficiently 
small,  only  one  polarization  mode  can  be  guided  by  the  fiber.  Then  we  have  single¬ 
mode  single-polarization  fibers. 

That  the  polarization  modes  in  the  conventional  single-mode  fibers  are  nearly 
degenerate  can  be  traced  to  the  fact  that  these  fibers  have  :  (i)  nominally  circular  cross 
sections,  (ii)  isotropic  index  and  (iii)  azimuthally  independent  index  profile.  When  one 
or  more  of  these  conditions  is  removed,  the  polarization  mode  becomes  non-degenerate. 
The  simplest  and  most  obvious  examples  are  fibers  with  noncircular  cross  sections  (and 
isotropic  and  (/(-independent  index  profiles).  Fibers  with  either  an  elliptical  core,  clad¬ 
ding  or  both  arc  of  these  category  and  have  been  studied  extensively  by  many  authors 
(3-6]. 

Fibers  can  also  be  made  nondegenerate  by  making  the  core,  cladding  or  both 
anisotropic.  If  either  region  is  biaxial,  a  fiber  becomes  nondegenerate  even  if  it  has  a 
perfectly  circular  cross  section.  Throtigh  mechanical  stress,  isotropic  materials  can  be 
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mncle  anisotropic.  The  stress  can  be  built  into  the  fibers  by  the  fabrication  processes  or 
applied  after  fabrication.  Of  particular  interest  is  index  anisotropy  induced  by  lateral 
stress.  Although  fibers  with  uniaxial  index  of  refraction  have  been  the  subject  of 
several  studies  [7-9],  fibers  with  biaxial  material  have  not  been  examined  until  recently 
(10-21].  Despite  of  these  studies,  fibers  under  lateral  stress  are  not  very  well  under¬ 
stood.  This  is  the  subject  of  this  work.  Recently,  a  series  of  papers  on  high 
birefringence  fibers  have  been  published  by  Snyder  and  his  associates  (16-21].  However, 
there  is  a  basic  difference  between  (heir  work  and  the  work  reported  here.  In  their 
work,  they  have  assumed  from  the  outset  that  LP|‘o  mode  is  coupled  with  LP^|  mode, 
and  LPfo  mode  with  LP^i  mode.  Here,  no  assumption  has  been  made  relative  to  the 
coupling  between  any  polarization  modes.  All  modes  are  included  in  the  consideration. 
The  inclusion  or  exclusion  of  a  particular  mode  is  determined  by  the  differential  equa¬ 
tion,  the  boundary  conditions  and  the  symmetry  of  the  problem,  and  does  not  rely  on 
our  prior  knowledge  of  the  field  distribution  in  the  fiber. 

In  the  following,  we  begin  by  reviewing  the  effects  of  the  stress  on  the  refractive 
index.  Under  lateral  stress, an  isotropic  medium  becomes  biaxial  (22,23).  Therefore 
waves  guided  by  biaxial  media  are  also  reviewed.  Complication  associated  with  biaxial 
media  are  noted  and  ways  to  solve  the  problem  are  described.  This  is  done  in  Sections 
2  and  3.  We  then  proceed  to  study  the  fields  and  dispersion  of  fibers  with  biaxial  core 
and  cladding  (Section  -1).  The  numerical  results  are  presented  in  Section  5.  Also  dis¬ 
cussed  are  the  similarities  and  differences  between  elliptical  fibers  and  biaxial  fibers. 
The  coupling  of  biaxial  fil)ers  is  discussed  in  the  last  section.  The  differences  between 
low-  and  liigh-  birefringent  fibers,  so  far  as  coupling  is  concerned,  is  also  noted  there. 

2.  Waves  Guided  by  Uniaxial  and  Biaxial  Media 

When  a  solid  is  stressed,  its  refractive  index  changes  and  this  is  known  as  the 
photof'lasl ic  elTecl.  Consider  an  isotropic  material  with  an  index  n  under  stress-free 
condition.  When  stressed,  its  indices  along  principal  stress  directions  become: 
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n,  n  +  [C,  T,,  +  C2(Tyy  +  T„)l 

n^n  +  {C,Ty,  +  C2(T,,  +T„)| 

^  n  +  [C,  T„  +  C2(T,,  +  Ty,)) 

where  Tyy,  and  are  the  stress  components  in  Cartesian  coordinates  and  C| 
and  C2  are  the  stress-optic  constants.  These  stress-optic  constants  are  negative  and 
are  of  the  order  of  10’’“  in  MKS  system.  More  importantly,  |  C2I  is  a  few  times 

larger  than  |C,|  (22,23). 

Since  the  fiber  core  and  cladding  regions  are  made  of  different  material,  the  ther¬ 
mal  expansion  coefficients  for  two  regions  are  likely  different.  After  fibers  are  drawn 
from  the  preforms  and  allowed  to  cool  down,  lateral  stress  is  frozen  in  the  fibers.  This 
is  the  basis  of  built-in  lateral  stress  in  many  polarization-holding  fibers  [24,25  and  26). 
If  in  addition,  the  fiber  cross  section  is  not  circular  and  we  have  /  Tyy.  In  the 

absence  of  any  build-in  axial  stress  (Tjj  =  0),  the  above  equations  become: 

“  n  +  (C,  T„  +  C,T„  I  ^  n,  -  (C,-C,)T„ 

",  “  n  +|C,  T„  +  CjT,.!  n,  -  (C,-C,)T„ 
n,  +1C2(T„  +  T„|| 

Clearly,  the  core  and  cladding  regions  are  bi<axial  with  n,  ^  n.  ^  n„. 

To  study  the  waves  guided  by  fibers  with  biaxial  core  and/or  cladding,  it  would  be 
necessary  to  examine  (he  waves  guided  by  biaxial  media  [27,28].  For  this  purpose,  we 
begin  with  the  time-harmonic  htaxwell’s  equations  for  biaxial  media  : 


V  X  E(x.y,z) 

=  ju,7<o  fl(x.y.z) 

(1) 

V  X  l!(x,y,z) 

=  f^r  ■E(x,y,z) 

(2) 

% 

'sc'i 


VH(x,y,z)=0  (3) 

V  c„rB(x,y,z)  =0.  (4) 

A  time-harmonic  variation  of  the  form  is  understood.  In  a  matrix  form  therelative 
dielectric  tensor  T  may  be  written  as 

n^  0  0 

r=  0  n^  0  (5) 

0  0  n2 

When  the  geometry  of  the  guiding  structure  and  the  relative  dielectric  tensor  are 
independent  of  z  and  waves  propagating  along  z  direction  are  of  interest,  a  term 
can  be  factored  out  explicitly.  Here  0  is  the  propagation  constant  of  the  guided  mode. 
It  is  also  convenient  to  separate  the  field  components  and  the  V  operator  into 
transverse  and  longitudinal  parts: 

E(x,y.z)  =  (l^Jx.y)  +  a3  e^fx.y)  (6a) 

ff(x.y.z)  =  [  K((x,y)  +  a,,  h^fx.y)  je*''*  (6b) 

V  =  V.  +  =  |i.  ■£-  +  i, +  j/?  i,  (7) 

where  is  a  unit  vector  in  +z  direction  and  the  subscript  t  signifies  the  components 

transverse  to  the  z  direction.  In  terms  of  these  variables  and  V^,  (1)  and  (2)  become 

V,  X  (8) 


K  X  (j'^t  -  ’^tCz)  = 


a,  X  (j^ht  -  Vth3)  =  -jw 


(11) 


where 


K  0 


By  eliminating  or  from  (9)  and  (11),  we  obtain 
=  j/^K-V^e^  -  ja;//„K-a3  x  V^h, 


Kt  =  jcj  a,  X  K  •  ft  •  -  j/?  a^  x  K  •  a,  x  h. 


whore 


k;2  0 

0  Ky2 


K  2  =  k2  n2  -  0^ 
K;  =  ko  n;  -  0^ 


and  k'o  -  oj"  fi^  fg.  DifTercntlal  equations  for  e,  and  can  be  obtained  by  substituting 
(13)  -  (H)  into  (3)  and  (4); 


K;-^  +1'; +  ''."'S-  >■.  =  -  ">)  (in 

Once  e,  and  are  known,  (13)  and  (14)  can  be  used  to  evaluate  the  transverse  field 
components. 

For  isotropic  as  well  as  uniaxial  media  (n,  =  ny),  clearly  K,  =  Ky  and  e,  and  h, 
arc  uncoupled  and  (10)  and  (17)  become  the  usual  wave  equations  for  e,  and  h,.  In  this 


-  6- 


fashion  ,  the  dispersion  relation  for  uniaxial  fibers  can  be  and  has  been  derived  (7,8,9). 
However,  the  difTerenccs  between  isotropic  and  uniaxial  fibers  are  noteworthy.  For  iso¬ 
tropic  fibers,  the  transverse  electric  fields  of  LPoj  mode  is  mainly  along  one  direction,  x 
or  y,  and  j'^tl  /|  h^|  is  approximately  a  constant.  In  uniaxial  fibers,  on  the  other  hand, 
both  e,;  and  Cy  are  present.  Each  components  of  and  have  to  be  expressed  as  the 
sum  of  two  parts,  one  involves  the  ordinary  index  of  refraction,  n,  =  ny  =  n,,  and  the 
other  the  extraordinary  index  of  refraction  n,  =  n^.  As  a  result,  |l?t|/|hj|  is  not  a 
constant. 

For  biaxial  media,  /  ny  /  n^^,  and  are  coupled  in  general.  Two- 
dimensional  waveguide  structures,  integrated  planar  waveguides  for  example,  are  excep¬ 
tions.  For  the  infinitely  large,  two  dimensional  structures  we  can  further  assume 
d  ^ 

=  0  or  =  0,  and  (16)  and  (17)  again  can  be  reduced  to  the  usual  uncoupled 

V 

wave  equations  for  and  hj.  This  is  precisely  the  approach  adapted  by  Burns  and 
Warner  in  their  study  of  uniaxial  planar  waveguide  structures  (24).  No  such 
simplification  is  possible  for  the  fibers  with  biaxial  media  and  herein  lies  the  crux  of  the 
problem.  In  addition,  and  Ky  can  be  positive  or  negative  depending  on  the  values 
of  kj,  n,(  and  k„  n^.  Thus  (16)  and  (17)  can  change  from  differential  equations  of 
elliptical  type  to  that  of  parabolic  type  as  0  changes.  Thus  for  biaxial  media,  (16)  and 
(17)  are  inhf’rently  difficult  to  solve,  and  no  exact  solution  is  known. 

Alternatively,  e,,,  and  h,  can  be  expressed  in  terms  of From  (4)  and  (6a)  we 


e*  -  Vt  (ri  •  (18) 

By  performittg  vector  product  of  on  both  sides  of  (8),  we  solve  for  a,  x  Simi¬ 

larly,  (0)  is  used  to  solve  for  a,  x  hi.  Upon  substituting  the  expressions  so  obtained 
into  (II)  a  difTcrential  equation  tor'll  is  obtained: 


o 


et  +  (k  2  n-  -  =  V J( V,  .  (r,  • 


Siniilnrly,  an  cli[Tcrontial  equation  for  is  also  obtained: 

+  (ko^<'j  +  0^  *2  X  V'  •  »2  X  )Kt 

=  ^i(^t  X  K^)  -  j  /?  a,  X  Vt"'  •  a,  x  V^h, 
In  component  form  (IQa)  can  be  written  .os: 


where  N  =  i3/kg  is  the  effective  index  of  refraction  of  the  guided  mode. 

AKhoiigh  and  are  still  coupled  in  (20a)  and  {20b),  the  differential  equations 
are  simplified.  In  particular,  for  the  fibers  of  interest,  |  n^  -  n^  |  /n^  <<1,  the  terms 
in  the  right  hand  side  of  (20a)  and  (20b)  arc  small  and  can  be  treated  as  perturbation 
terms.  Then  approximate  solutions  can  be  obtained  from  (20a)  and  (20b).  This  will  be 
done  in  .Section  i.  However,  some  features  of  the  solution  can  be  deduced  from  (20a) 
and  (201))  without  actually  solving  the  differential  equations.  Clearly  if  e,  is  an  even  (or 
odd]  function  of  x  or  y,  then  Cy  must  be  odd  [or  even]  function  of  x  or  y.  When 
expre.ssed  in  tc'rms  of  cylindrical  coordinates  r,<f>,  e^,  and  Cy  can  be  written  as  Fourier 
sine  and  cosine  series.  Thus  for  the  core  region  (r  <  a), 


+  EK'i'f'')  ^x’i''('')  sin  i<^| 


ey(r,<?:))  -  eJ])(r)  +  V[e!|''(r)  cos  i</>  +  e^'(r)  sin  i.i!>| 


(21a) 


(21b) 
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and  for  the  cladding  region  (r  >  a): 

ex(r,^)  =  exo(r)  + 

i  =  l 

cos  \4>  +  ejJ]''(r)  sin  i^] 

{22a) 

_  00  _ 

ey(r,^)  =  ey^o(r)  + 

i=l 

cos  \<j)  +  e^^'(r)  sin  \<f>] 

(22b) 

The  subscripts  R  and  D  are  used  to  signify  that  these  terms  are  associated  with  the 
core  and  cladding  regions  respectively.  From  (20a)  and  (20b)  it  is  clear  that  the  sine 
and  cosine  series  of  are  coupled  respectively  with  the  cosine  and  sine  scries  of  Cy.  In 
addition,  the  even  harmonics  are  coupled  only  with  even  harmonics,  and  odd  harmonics 
with  odd  harmonics  only. 

3.  Fields  and  Boundary  Conditions  for  Weakly-guiding  Fibers 

An  isotropic  fiber,  with  core  and  chadding  indices  n,  and  nj  respectively,  is  viewed 
as  "weakly  guiding"  if  n^  -  nj  «  n,.  Under  the  condition  of  weakly  guiding, 
and  I  V,hj|  ,  are  negligibly  small  in  comparison  with  |  and  |  hj]  ,  and  therefore 
can  be  neglected.  Then  from  (9)  and  (1!)  we  have: 

hi  -  ajiX^i  (23) 

which  is  accurate  to  the  order  of  |  (n^  -  nj)/n,,|  .  Thus  for  weakly-guiding  fibers, 
and  h|  are  related  in  the  same  manner  as  that  of  uniform  plane  waves  in  isotropic 
media  [29].  This  is  a  basic  feature  of  fields  of  weakly  guiding  isotropic  fiber. 

Now  consider  weakly  guiding  anisotropic  fibers.  In  the  core  region,  the  indices 
along  three  principal  directions  are  n^,,  n^y  and  n^j.  Correspondingly,  the  indices  in  the 
cladding  region  are  n^„,  n^y  and  n^^.  The  core  radius  is  a.  By  weakly  guiding,  it  is 
meant  that  |  (n,,^  -  n,|,()/nr^|  etc.  are  small  and  are  of  the  order  of  0(f).  Under  these 
conditions,  following  the  same  argument  presented  by  Gloge  for  isotropic  fibers  [29],  we 
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(24) 

With  known,  oxaci  expr<ssions  for  and  e,  and  approximate  expression  for 
can  be  obtained  from  (9).  (18)  and  (21). 

It  has  been  noted  in  the  last  section  that  the  cosine  series  [and  sine  series]  of  Cj^ 
and  sine  series  [and  cosine  series]  of  Cy  are  coupled  through  the  dilTerential  equations 
(20a)  and  (20b).  We  like  to  show  that  the  same  conclusion  can  also  be  reached  by  con¬ 
sidering  the  boundary  conditions.  To  enforce  the  boundary  conditions,  it  is  expedient 
tf)  express  and  in  terms  as  e^,  e^,  and  In  terms  of  the  coellicients  intro¬ 
duced  in  (21a)  and  (21b)  the  field  components  in  the  core  region  are: 

2e^  =  -  e]i'=) 

+  (e|?/  -  2efo  +  ey'l")  sind  +  (ey^2^  +  2ey^  -  0^2!  cosfii 
£  h^xYi  +  +  Cy'^j-l  +  Py^-  +  i)  sinj<A 


SV 

+  S  I^X^'-I  -  COSj(^ 

j-^ 


=  4t-'.V'  +  T  <-.Vl  +  ^  c«,-| 


n"  t  n  r' 


'I?  '  '-■  ’T.F  l'J,P''C"’Vr  ’’-• 
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2/? 

-  -  ' 

wf„ 


O  f?  j*  I 

-  hr“xPxl  +  n 


,2pRsl 

‘ry^y  1 


+  Inrxex’2'  +  2n2e]^o  +  ^ry^'y'll  ‘'0S<A  +  [n«ex^2®  +  2n ~  nj Cy’^al  sin(^ 
+  i;[nr^xCx^-i  +  n^c^f+j  -  njej^|_,  +  njej^?+,  jcosj^ 


+  Slnr^xex’^J-i  +  n^ej^l+i  +  nJefjL,  -  nJeJ^f+j  |sinj(^  (28) 

j=2 


In  these  expressions  a  prime  indicates  d/dr.  Corresponding  expressions  can  be  written 
for  the  field  components  in  the  cladding  region,  with  the  superscript  R  replaced  by  D 
and  n^x  etc.  by  nj,<  etc. 

At  the  core-cladding  interface  (r  =  a),  e^,  e^,  d„  h^,  h,  and  b,  should  be  con- 
tinous.  It  is  known  that  when  the  continuity  for  the  tangential  components  of  fi  and 
ft  are  satisfied,  the  continuity  for  the  normal  components  for  0  and  ft  are  automati¬ 
cally  satisfied  [30]  Since  sin  ](1>  and  cos  j<f)  are  orthogonal  functions  in  the  range  (0,  2a-). 
the  continuation  of  the  tangential  components  of  ft  and  ft  amounts  to  the  continuation 
of  the  individual  Fourier  components.  More  specifically, 


(29a) 


2  2 

+  1  „Rcj  +  illLfpRs/  .  J.  -Rs] 

o  l‘xl  ^  „  ^xll  ^  2  ^  ' 


+  -  o.'i-i  +  -  '"’I 

"Jz  ^  "dz  ^ 


(29b) 


Rr  1  ^  f?«;  —  *>  Dr  I  *> 

+  flryVl  '  "dx^  I  +  «dy^: 


cly'  y  I 


(29c) 


^x2  “^xO  ^  ^y2  ^x2  ‘‘'^xO  ”y2 


,Rc  4.  On^  —  4-  9o^  — 

y2  '  ^^yO  ^x2  ^y2  '  ^^yO  ^x2 


■^(2<-.V  +*'  +f’-x"/)  +  -f '<■»"’’  +f  ef2’) 


■(2cft'  +<>?2''  +-4')  +  ^(<'"/'  +7®,'^) 

a  n|,  a 


2  2 
i!£L/„Rs»  4.  A„Rs»  4.  il£Lf9„R»_„Rc#_ApRc\ 
2  'f^x2  +  „  <^x2 )  ^  2  y2  o  y2  > 

n„  a  n„  » 


njz  ®  "dx  ^ 


n2„Rc  4_  nn2pR  4.  n2oR<  —  n2,,Dr  4.  9  2pD  4.  _2pDs 
^rx^x2  ^  "^rx^xO  ^  ^r>'^y2  **dx'x2  ^  ‘'*^dx®xO  ^  ^dy®y2 

„2pRs  4.  9,,  2pR  _  2„Rr  -  n2  pOs  4.  9n2  pR  _  n2  pPc 
Rrx^x2  ^  *Rry^yO  Rry^y2  *^dx^x2  “^dy^yO  Rdy^y2 

a(o,'!,''  -  2..;|,'  +  oK,"')  +  2(cS'  +  0,1') 

=  2(0.1'  -  2'-.'o'  +  '•y'o" )  +  2(4'  +  4') 

a(e1'  +  2o,V  -  0,1' )  +  2(0,1  -  o.1l 

=  2(0,1'  +  2o,1  -  0.1' )  +  2(0,1  -  0,1) 


And  for  j>2, 


I  -ILuRr.  I  4..Rc  I 


^  o  (®yj+l 


2 

-  f  ^/„Dc  I  I  „Dc  f  +  i±lpPc  _jlJLpPr  \ 

-I  o  l^xj-l  ^e^j+1  ^  ,  <^xj  +  l  „  e^j-l/ 

T\  **  a,  A 


2 

+  fpRs  /  _pPs »  j.  jzi-pP?  a.  J-iLLpP?  u 


f  illlrpRs  I  4-pR-''  / -U-p R.'^  .  ■+  J-tJLp Rg, . ] 

I  o  1  ^^xj  +  l  „  ^xj  +  l^  „  6xj  +  |) 

n  *•  di  a. 

"rz 


2 

4-  ^^y  fpRc  f_oR«  >  J~^  I, Rc  J  ^  p Rc  \1 

^  2 '^yi-1  ^yi  +  i  o  ®yri  o  ®yi  +  »'J 


!liiL,„Ds  I  4._D!.  *_I_LpPs  1_  j  +  1  Ps  \ 
o  l^xi-1  ^^xj  +  1  _  ^xj+l^  _  ®xj  +  l/ 


n  l^xj-1  ^xj  +  1  ^xj+l^  ^xj  + 

njz  ^  ^ 


4-  ili^fpRc  f_pDc  f_JzLpDc  - j  tJ-pPc  \1 

“dz 


n2„Rr  4.  „2_Rr  _  p2pRs  4.  2pRs 

Rrx*^x  j  I  ^  j  +  1  ^ry^y  j-I  ^ry^y  j  +  I 


—  n2„Dc  4.n2„Dc  _n2-Ds  .|._2_Ds 
^dx^x  j-|  ^  *^dx®xj  +  l  ^dy*^y]-l  '  ^^dy'^yj  +  l 


(31c) 


(3  Id) 


j-1  "rx^xj  +  I  “rv'yj-l  ”ry^yj  +  i 


=  "d^x^j-i  +  "dx^i^j  +  i  +  -  ndVy'’j  +  i 

+  ei'iV.'  +  eyY.'  +  ^y^i  +  i'l 

+  [(rl)Cx’^j-i  +  (j  +  Oex’^[+i -(rl)ey’^^l  +  (j  +  l)ey’^f+i] 

=  a[-eP[./  +  cDf_/  +Cyf>f„'] 

+  ((rO<^x^'-i  +  (j  +  iK^^i  -  (i-i)<^y'^j-i  +  (j  +  i)ey’^j+il 

a[eJ^^l'  -e]^f+l'  +^yi-i 


(31g) 


+  (  -  (rl)Cx^'-i  “  (i  +  OCx'^j  +  i  -  (j-l)ey'^f-i  +{j  +  l)ey’^/+i 


^l^x  j-1  ®x  j  +  t  ^  ®y  j-l  ®y  j  +  1  i 


+  [  -  (j-llei^j-i  “  (j  +  lK'^j  +  i  “  (j-l)ey'^[-i  +(j  +  l)ey‘^f+il 


Ds  _  .  _  li_iuDc 


(31h) 


rAamination  of  these  boundary  condition  equations  reveals  that  the  cosine  series  of 
Gj,  are  coupled  with  the  sine  series  of  e^,  and  the  sine  scries  of  e,,  with  the  cosine  series 
of  Gy.  Furthermore,  the  even  harmonics  are  coupled  only  with  the  even  harmonics,  and 
odd  harmonics  with  odd  harmonics.  For  example,  in  (30a)-(30d),  and  e;jo  are  cou¬ 
pled  with  e,[2'’  ^y*2’'  ^x^2'  ^y*2^  through  these  terms  to  the  higher  order  even  har¬ 

monics. 

In  the  ca.se  of  isotropic  fibers,  ej[o  and  e^  |or  and  e^]  are  needed  to  determine 
the  fields  and  dispersion  of  LI’oi  mode,  No  higher  order  terms  is  required.  It  would  be 
interesting  to  see  if  the  same  is  also  true  for  biaxial  fibers.  Suppose  that  only  one  term 
from  each  series  is  kept,  (i.e.,  e^^  for  the  core  region  and  ej^  for  the  cladding  region). 
Then,  the  boundary  conditions  (30a),  (30c), (30e)  and  (30g)  become 


®x0  ®x0 

(30al) 

r.2  n2 

"rx  R'  _  "dx  D' 

2  ^xO  9  ®x0 

(30c 1) 

^rz  ^dz 

_2  pR  -  ^2  pD 

^rx  ®x0  "^dx  °x0 

(30el) 

o'  =  e'  ^9 

^  xO  ^  xO 

(30gl) 

For  isotropic  fibers,  n^,,  =  n„  =  and  n^^  =  nj„  =  n^,  (30cl) 

is  identical  to  (30gl). 

Rccnll  the  expression  for  is  exact  while  that  for  is  approximate  and  accurate  to 
0(c).  When  terms  of  the  order  of  0(c)  or  smaller  are  ignored,  (30el)  also  reduces  to 
(30al).  Therefore  in  (he  case  of  weakly  guided  isotropic  fibers,  the  boundary  conditions 
are  simply  the  continuation  of  c„o  and  e\o 

For  biaxial  fibers,  n^^  ^  n„  (30cl)  and  (30gl)  are  not  identical.  The  pres¬ 
ence  of  etc.  makes  it  impossible  to  satisfy  all  boundary  conditions  if  and 

are  the  only  terms  kept.  There  have  to  be  the  higher  harmonic  terms. 

Suppose  only  the  leading  harmonic  term  of  cosine  series  of  c,  (e,*Q  and  e^)  and  the 
sine  scries  of  Cy  (Cy*!"  and  Cy^'’)  are  kept.  There  are  two  unknown  Fourier  coefficients  for 
the  fields  in  the  core  region  and  two  in  the  cladding  region.  In  addition  the  propaga¬ 
tion  constant  0  is  also  an  unknown.  Four  equations  [(30a),  (30c),  (30e)  and  (30g)] 
are  used  to  eliminate  four  unknowns.  The  result  is  a  dispersion  relation  determining  P 
or  the  offeciive  index  N.  If  desired,  we  can  use  the  first  two  or  first  three  harmonic 
terms  of  the  sine  and  cosine  series,  then  we  have,  respectively  8  or  12  unknown 
coefficients.  Thus  8  or  12  boundary  condition  equations  respectively  arc  needed  to  set 
up  the  dispersion  relation  for  /?.  Dy  solving  the  dispersion  relation,  we  can  determine  /} 
as  a  function  of  various  fiber  parameters.  In  addition  all  but  one  unknown  Fourier 
coefficienfs  can  be  expressed  in  terms  of  one  Fourier  coefficient.  The  undetermined 
Fourier  coefficient  can  be  determined  from  the  field  strength  of  the  mode  or  the  power 
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carriod  by  the  fields. 


4.  Approximate  Expressions  of  e^  and  e^  and  Dispersion  Relation 

To  derive  the  approximate  expressions  for  the  transverse  electric  field  components, 
we  treat  the  right-hand  terms  of  (20a)  and  (20b),  as  perturbation  terms.  Thus,  by 
ignoring  the  terms  on  right-hand  side  of  (20a)  and  (20b),  we  obtain  the  zeroth-order 
solutions,  ej°*  and  With  the  zero-th  order  solutions  substituted  in  the  right-hand 
side  terms  of  (20a)  and  (20b)  and  treated  as  known  quantities,  then  (20a)  and  (20b) 
become  inhomogeneous  differential  equations.  The  particular  solutions  for  these 
differential  equations  are  then  combined  with  the  homogeneous  solutions  which  are 
exactly  ej®'  and  Thus  we  have  and  which  are  accurate  to  0(f). 

Two  lowest  order  polarization  modes  are  mutually  orthogonal  to  each  other.  In 
one  polarization  mode,  its  e,,  is  described  by  a  cosine  series  and  Cy  by  a  sine  series.  For 
the  other  polarization  mode,  and  Cy  are  given  by  sine  and  cosine  series  respe(  vely. 
In  this  section,  our  attention  is  focused  on  one  of  the  polarization  modes,  namely  the 
mode  with  a  cosine  series  for  c„  and  a  sine  series  for  Cy.  Since  no  restriction  is  placed 
on  the  value  of  n,„  relative  to  n^y,  or  nj„  relative  to  njy.  The  other  polarization  mode, 
can  be  obtained  simply  by  interchanging  n„  and  nj,  with  n^y  and  njy. 


Fields  in  the  core  region; 

To  calculate  the  fields  in  the  core  region,  we  introduce  two  new  variables, 
x"  =  (n^y/n^Jx,  and  y'  =  (n„/n„)y.  Thus  (20a)  and  (20b)  become: 


+  -tK'  (n  “  -  N"^)  e,  = 


n„n„  dxdy' 


(32a) 


When  the  right-hand  side  terms  in  (32a)  and  (32b)  are  ignored,  we  have  immediately 
the  zero-order  solutions: 


n  _>  00  n,, 

eio)  =  Ao  .Io(—  U,,  ^)  +  EA.  U-f-  Urx  t) 


n„  ■■■  a 


"r*  a 


00  n  ^ 

ef'  =  I;B„  J„(-^  sin  n,^” 


■ry 


where 


U„  =  k„a  \/n2  -  n2 


Ury  =  k„a  y/ n  J  -  N“ 


r'  =  (x2  y'2]-</2  =  (x2  +  ^  <!>'  =  tan->(-5^^) 


Dr,X 


rM 


=  [x"2  +  y2]  =  +  y^l  (i>"  =  tan'’(-^) 


*r2 


nryX 


(33) 


(34) 


(35a) 

(35b) 


(35c) 


(35d) 


As  stated  previously,  only  cosine  and  sine  scries  are  kept  in  (33)  and  (34)  respec¬ 
tively.  In  defining  the  new  variables  y'  or  x" ,  the  coordinates  are  scaled  by  constant 
factors  but  there  is  no  translation  in  the  origin.  Thus  the  symmetry  properties  of 
and  Cy  with  respect  to  \" ,  y' ,  or  remain  the  same  as  those  with  respect  to  x 
and  y,  and  (f),  as  discussed  in  Section  3. 

With  e]*’*  and  e]®'  known  and  substituted  into  (32a)  and  (32b)  respectively,  (32a) 
and  (321))  become  a  set  of  inhomogeneous  differential  equations.  Their  solutions  are 
comprised  of  two  parts:  a  homogeneous  solution  and  a  particular  solution.  The  homo¬ 
geneous  solutions  are  of  the  same  form  as  the  zeroth-ordcr  solution.  In  view  of  the 


7^ 

IV-.i 

’••.■•I 


k- 
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functional  form  of  and  ej®^,  we  can  expect  that  the  particular  solutions  can  also  be 

Drj  f  r 

expressed  as  series  of  J^( - LT  — )  cos  m  ^  or  J^{ - —)  sin  m  <f>.  By  substitut- 

Dry  a  n„  a 

ing  these  series  into  the  left-hand  side  of  (32a)  and  (32b),  performing  the  necessary 
differentiation,  and  comparing  the  resulting  terms  with  those  on  the  right-side,  the  par¬ 
ticular  solutions  are  found.  In  the  lengthy  manipulation,  the  recurrence  relations  for 
Bessel  functions  [31]  are  used  repeatly  and  we  also  note  in  particular  that  U„  U^. 
Then  ej**  and  ej’\  accurate  to  0((),  arc  obtained  by  combining  the  the  particular  solu¬ 
tions  with  the  zcro-th  order  terms: 

ei‘>  =  A,  Jo(—  Urx  -)  +  EA„  -t) 

n,,  a  ,  n„  a 

+  1  "ry(nr^rn^) 

4  nrz("r^x“nj) 


E  ^)  COS  (m-2)^  -  Jm.f2(-~  U,y  COS  (m  +  2)(^  ] 


•ry 


*ry 


!"=  K,,-!— )sinmf 


ft 


*ry 


1  n„(n.Vn.^,)  n„  , 

+  -  V-  '  -  ,  {  2A,  J,!—  U,,  sin  U 

4  n,.(n2-„2|  n4  '  ^  *'n„  "  a' 


(36) 


+  £a„  U„  ^)sin  (m+2)(i  -  U„  -f)  sin  (m-2)^|} 


1  *^rx  a  Hf,; 


(37) 


Note  that  |  n  )/[n„(n2-n2j)|  and  |  n^yfn^-nj )/(n„(n2-n2)]|  are  of  the 

order  of  0(1),  while  |  (n,^jj-N^)/N^|  and  |  (n^-N^)/N^|  are  of  the  order  of  0((). 


Fields  in  the  cladding  region: 

Again,  new  variables  may  be  introduced  to  simplify  the  differential  equation.  In 
particular  r' ,  ,  r' '  and  <{>' '  are  defined  in  the  same  fashion  as  given  in  (35c)  and  (35d) 

with  n,,(  etc.  substituted  by  n^,;  etc.  For  isotropic  fibers,  N  of  a  guided  mode  is  always 
larger  than  nj.  For  biaxial  fibers,  complication  arises  particularly  for  modes  near  cutoff 
since  n^,;  7^  njy.  Near  cutoff,  it  is  possible  to  have  nj„  <  N  <  njy  or  nj,  >  N  >  njy 
depending  on  the  relative  values  of  and  njy.  Thus  three  cases  have  to  be  treated 
separately.  For  convenience,  we  define: 


Wdx  =  k,a  v/N^  - 


(38a) 


W^y  =  k,a  v^N^  -  n|y 


(38b) 


U,.  =  k„a  x/nl  -  N2 


lJ,y  =  k„a  \/nly  -  N- 


(38d) 


For  N  >  n(j„  and  nj^,  and  Wjy  are  positive.  Following  the  same  procedure 
described  in  previous  section,  we  obtain: 

e'l)  =  C„  K„|-^  Wj.  +  £c„  K„(^  Wj,  )  cos  m,S' 

^dx  ^  I  '’dx  a 


1  ndv(ndVn,iv) 
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g  D„  I  Wj,  i)  eos  (m-2)0  -  Wj,  f)  cos  (m+2)^  ] 

T  ^dy  ^  ^dy  a 


/  .  00  rij-  r' ' 

el')  =  D^n,  K„(^  Wj,  J— )  sin  m^' 

I  "dy  ^ 


1  ^dx(^dz  ^d"^x)  ^dx  ^  ,,  ,  "^dx  ,,,  ^  \  •  />  ; 

7  TTT-nr  “7^^  2Co  M—  w,,  -)  sm  2<A 


+  EC„,  1K,+2(-^  W,,  sin  (m+2).^  -  K,.2(-^  W^,  sin  (m-2)<il} 

1  “  dx  ^  ”  dx  ^ 


If  n^x  >  n^iy,  there  exists  situations  where  nj^  >  N  >  njy.  With  >  0  and 

Ujy  >  0,  the  fields  in  the  cladding  are  described  by  a  mixture  of  the  Hankel  function 
and  the  modified  Itessel  function.  The  presence  of  the  Hankel  function  signifies  the 
exi.stence  of  the  leaky  waves  in  the  cladding  region. 

e!"  =  f'o  I'ol-^  Wd.  -t)  +  gc„  K„(-^  W,,  .^)  cos  m«' 

"dx  a  ,  nj,  a 

j_  ndy(»,Tz-»dy)  n,yy-N- 

^  nj,(njVndy)  N‘ 


’2(“^  ^’dy  ~)  (m-2)<;i  -  Hdy  “)  (m  +  2)(>  ) 

ndy  a  rijy  a 


(41) 


1_ 

■*  "d.(oi-“d^)  N2 


{  2Co 

“dx 


Wdx 


— )  sin  2<f> 
a 


oo  I1j_  r  Dj-  p 

+  EC„  1K„+2(—  W,,,  i)  sin  (m+2)^  -  K„-2(—  Wj,  i)  sin  (in-2W) 
1  njx  a  n^x  a 


(42) 


On  the  other  hand,  if  n^^^  <  n^^,  there  may  be  situations  where  U|,j  >  0  and 
>  0.  Then  the  fields  are  given  by: 


,l'l  =  C„  \m—  Ud.  -)  +  sc^  Ud.  )  cos  m^' 

"dx  “dx  ^ 


I  ^dy)  ^dy  ^ 

4  Odxind^-ndV 


£  u„  (  K„.2(—  Wd,  -)  cos  (m-2)4>  -  K„+2(-^  Wj,  f)  cos  (m+2)«  | 

1  I'dy  ^  “dy  ^ 

(43) 

oo  n  j,  r' ' 

ej')  =  KJ-^  Wd,  — )  sin  m^" 

1  "dy  ^ 


ndx(nJ2-nirx)  nj,-N^ 

»d.('>d,"”d^i 


{  2C„ 


^dz 

Ddx 


Udx  — )  sin  2<j> 
a 


+  I",a2(—  Ud,  -)  sin  (m+2)«  -  Uj,  f )  sin  (in-2)0|l 


(44) 


r  '  • 


IS.'/ 

4V  ,> 


•i 


Further  approximations: 

Before  imposing  the  boundary  conditions,  it  is  necessary  to  revert  r' ,  r"  ,0'  and 
back  to  r  and  <l>.  For  points  in  the  core  region, 


r'  -  [x“  +  -^y“]^‘/‘  r  -  -7(1 — ^)(l-cos2<?) 


(45a) 


n2  n2 

•Irv  «  n.  •  r  U. 


r"  —  [~^x^  +  y")  lii  r — ^(1 — ^){l+cos2(^) 
n“  4  n“ 


(45b) 


0'  ~  0  — ^(1 - —)  sin20 

2  n„ 


(45c) 


and. 


0"  ~0 — ^(1 — ^)  sin20 

2  Dry 


(45d) 


Thus,  to  the  order  of  0(e), 


m  , . 


cos  m0'  cs;  cos  m0  +  —  (1 - )  |  cos  (m-2)0  -  cos  (m+2)0 

4  n„ 


{45e) 


m  “rj 


sin  m0"  ~  sin  ni0 - (1 - )  (sin  (m  +  2)0  -  sin  (m-2)0]  (45f) 

4  n^ 


Making  use  of  these  approximations  in  conjunction  with  Neumann’s  addition  theorem 
for  the  Bessel  functions  [31]  and  keeping  only  terms  of  the  order  of  I  and  e,  we  have 


!i 

a 


r  >_1  "fL  IT 


U—  ^rx  -)  -  Jm(—  ^rx  -(“T  ^(  1"^ )( I-COs20)[J' J  7^  U,,  ^)| 


a  4  n 


a  n  “ 

"rz 


h 

n. 


•  “>  -  -  ’  -  ■ 


^  ^  ».  ^**4  «  m.r  ^  J1  ,  .  ■  tkJ*.  tkJ'  JX./.  V  •*» 
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“  ■’"•'■ir  f '4  4  i(‘-4h‘ +™s2«)|j'  u„  t'l 

Ufy  d,  llfy  d,  ^  lljy  d  JJ  ^  ^  U  j  y  H 


r*  I  r  r 


where  J'  indicates  the  derivative  of  the  Bessel  functions  with  respect  to  its  argument. 
Use  these  approximations  and  sorting  out  the  terms  with  different  harmonics,  we 


ej"  “  A,.  I  J«(-^  +  U„  e„  i)) 

Hr,  a  4  n,,^  a  n,,<  a 


+  Ao 


1  ^rz  r  ^r*  F  2  ^r*  f 

i  'if  2(7^  7)  +  f)i 

8  n^,(  a  a  4  n^^  Uj.,!  a 


n«2  r  1  ^rz  r  ^rv  ^rz  r 

+  A.  (  U-f-  U„  7)  -  7  7^  u,,  7(1-— )j' 2(7^  u„  f )  I  cos2^ 


+  A2  1 


+  A,  [ 


^  8 

“rx 

rx  -( 

( 

a 

“rz 

/  ”  I  T 

{ 

'-'rx 

“rx 

1 

nrz  , 

8 

s 

X 

/• 

1 

"rz 

ot 

"rx 

^rz 

f  1  T 

l 

»rx 

'^'rx 

1 

!!ii, 

‘rx  “rx 


a  4  n. 


a  n; 


a  n; 


"rx  a 


r  ^rx  .  ^ra  T  2  ^r»  ^ra  T 

u„  -(1— ?-)J'j(—  u„  -)  -  7(1-— )Jj(—  u„  -)  I  C0S4^ 
a  n,,  a  4  n„  n„  a 


^rz  r  1  ^ra  F  ^rx  ^ra  F 

+  A,  [  J,(^  U,,  7)  -  7  7^  U,,  -(1— 7)J7(-^  U„  ]  cos4«i 


a  4  n. 


a  n; 


a  n; 


^rx  a 


Tv  4  .  f 

7^  ■  7<*“7~^'^4{7-  ^'rx  7 

a  4  n,,  n,,  a 
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I 

-1^6  (  T 


8  n 


ry 


ry 


+  - 


1  nrx(nn"nr‘x)  n«-N2 


4  n,,(n2-n2j 


(2A0-A4J  J2(-^^U„-^^)sin20 
n. 


+  A2  I  J,(-^  U„  i)  sin#  I 


+  A,  I  J.(-^  U„  i)sin#| 


+  ... 


(47) 


By  rearranging  the  terms  in  (46)  and  (47)  in  the  forms  of  (21a)  and  (21b),  we  have  the 
Fourier  coeiricients  for  the  core  region, 


=  Aof(,o(r)  +  A2fo2(r)  +  B2go2(r)  (48a) 

eJI"!'’)  =  Aofjoir)  +  A2f22(r)  +  A4f24(r)  +  B4g24(r)  (48b) 

<^x’4'('')  =  A2f42{r)  +  A4f44(r)  +  B2g42(r)  -b  B6g46(r)  (48c) 

=  B2P22(r)  +  B4P24(r)  +  A^qsoCr)  +  A4q24(r)  (48d) 

ey'^lr)  =  B2P42(r)  +  B4P44(r)  +  B6P4g(r)  +  A2q42(r)  (48e) 

Explicit  expressions  of  foo(r)  etc.  are  given  in  Appendix  A. 


Similarly,  from  (11)  -  (41),  we  have  for  the  cladding  region: 

~  ^0  ^oo(*’)  ^2  ho2(^)  ^2  ko2(r)  (49a) 

»^x^2'(r)  =  Co  h2o(r)  +  C2  h22(r)  +  C4  h24(r)  +  D4  k24(r)  (49b) 


f 


=  C2  h42(r)  +  C4  h44(r)  +  D2  k42(r)  +  Dg  k46(r) 


(49c) 


Cy‘i‘'(r)  =  D2  r22(r)  +  D4  r24(r)  +  Co  S2o(r)  +  C4  S24(r) 

ey‘4(r)  =  ^>2  '■42(0  +  1^4  r44(r)  +  ^>6  »‘46('‘)  +  Cg  S42(r) 
Expressions  for  hoo  etc.  are  also  listed  in  Appendix  A. 


(4gd) 

(49e) 


Dispersion  relation: 

Three  leading  harmonic  terms  of  the  even  cosine  series  of  e,  and  the  even  sine 
series  of  are  kept  in  (41),  42),  (50)  and  (51).  The  coefficients  A^,  A2,  A4,  B2,  B4  and 
Bg  are  pertaining  to  the  core  region  and  Cg,  C2,  C4,  D2,  D4  and  Dg  to  the  cladding 
region.  Substituting  these  equations  into  the  boundary  condition  equations  (32a),  (32c), 
(32e),  (32g)  and  (33a),  (33c),  (33e)  and  (33g)  with  j  =  3  and  5.  we  have  12  linear  equa¬ 
tions  which  can  be  cast  in  matrix  form; 


Qu 

Q.2  • 

■  QlO 

Qla 

Qib 

Qic 

^0 

Q21 

Q22  ■ 

■  Q2fl 

Q2a 

Q2b 

Q2c 

Ao 

Qai 

Q32  • 

■  Q39 

Q3a 

Qsb 

Qsc 

A4 

Q41 

Q42  ■ 

Q49 

Q4a 

Q4b 

Q4C 

B^ 

Qsi 

Q52  • 

Q59 

Qba 

Qsb 

Qsc 

B4 

Qbi 

Q62  • 

•  Qe9 

Qoa 

Qeb 

Qdc 

Bs 

Qti 

Q72  • 

■  Q79 

Q7a 

Q7b 

Qic 

Co 

Qsi 

Q82  • 

•  Q89 

Qsa 

Qsb 

Qsc 

C2 

Q91 

Q92  • 

•  Q99 

Q9a 

Q9b 

Q9c 

C4 

Qal 

Qa2  ■ 

Qa9 

Qaa 

Qab 

Qac 

1)2 

Qbl 

Qb2  • 

■  Qb9 

Qba 

Qbb 

Qbc 

1)4 

Qci 

Qc2  • 

■  Qc9 

Qca 

Qcb 

Qrc, 

Expressions  for  Qjj  with  i  or  j  =  1,  2,  3,...,  0,  10,a,  b,  c.  are  given  in  Appendix  B. 

A  dispersion  relation  is  obtained  by  setting  the  12  x  12  determinant  to  zero. 

II  Q  11=0  (51) 

Obviously,  (51)  will  be  solved  numerically. 
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5.  Discusssio  and  Numerical  Examples 

Before  numerical  results  are  presented,  it  would  be  helpful  to  enumerate  the  simi¬ 
larities  and  differences  between  elliptical  fibers  and  biaxial  fibers  purely  from  a 
mathematical  view  point.  As  discussed  in  Section  4,  (20a)  and  {20b)  can  be  reduced  to 
the  usual  wave  equations  by  a  simple  substitution:  x"  =  (n^^/n^)  x  or 
y'  =  (n„/n„)  y  for  the  core  region,  and  corresponding  substitution  for  the  cladding 
region.  Thus  a  biaxial  fiber  with  a  circular  cross  section  can  be  viewed  as  an  isotropic 
fiber  with  an  elliptical  cross  section  so  far  as  the  differential  equations  for  the  zero-order 
solution  are  concerned.  However  the  scaling  factors  for  the  core  and  the  cladding 
regions  are  not  necessarily  the  same,  so  the  ellipticity  of  the  core  may  differ  from  that 
of  the  cladding.  The  main  difference  between  two  types  of  fibers  is  the  way  e,  and  Cy 
are  coupled.  For  elliptical  fibers,  the  coupling  occurs  only  at  the  core-cladding  boun¬ 
dary.  For  biaxial  fibers,  coupling  originates  not  only  from  the  boundary,  but  also  from 

the  entire  biaxial  region  or  regions.  This  is  due  to  the  terms  and  •— — —  in  (20a) 

axoy  axdy 

and  (20b).  Even  for  the  coupling  at  the  boundary,  there  is  a  difference.  In  matching 
df  at  the  core/cladding  boundary,  it  is  necessary  to  account  for  the  difference  between 
n^f  and  nyp  nj,j  and  ny^. 

Since  the  characteristics  of  optical  waveguides  depend  more  on  the  index 
differential  and/or  index  anisotropy  than  the  index  itself,  the  results  can  best  be 
described  in  terms  of  index  differential  and  anisotropy.  As  stated  previously,  the 
indices  of  the  core  and  cladding  regions  are  (n„  ,  n^y  ,  n„)  and  (nj^j  ,  njy  ,  n^,)  respec¬ 
tively  (Figure  1).  The  index  difference  between  the  core  and  the  cladding  regions  is 
defined  as 

^rd  ~  ^^drnl/'^rm  (^2) 


where 


iifm  =  max  of  (a,,,n,y) 
njm  =  of  (Ujjj.Ildy) 

Similarly  the  index  anisotropy  of  the  core  and  the  cladding  regions  are: 

^xyr  —  (^rx  ~  ®ry)/^rx 


^xzr  (^rx  ®rx)/**r 


(53b) 


^xyd  (*'dx  ^dy)/^dx 


^xzd  (^dx  ^dzl/^^dx 


(53d) 


In  all  numerical  calculations,  nj^^  is  set  to  1.5000  while  other  indices  are  varied  by 
changing  d„yp  d^yj  and  d^j  etc.  To  make  the  numerical  results  more  meaningful,  they 
are  cast  in  terms  of  generalized  parameters  V  and  b, 


V  =  ka 

b  =  (N2  -  n^„,)  /  (n^^  -  n|^) 


(54a) 


(54b) 


which  are  the  same  as  those  used  for  isotropic  Dbers. 

The  transverse  electric  fields  for  two  lowest  polarization  modes  have  been 
obtained,  using  the  method  developed  in  the  previous  sections.  The  fields  for  the 
lowest  polarization  modes  are  depicted  in  Figures  2a  and  2b,  and  those  for  the  second 
lowest  polarization  mode  are  given  in  P'igures  3a  and  3b.  In  these  plots,  |e,(|  and  |  Cy] 
are  plotted  as  functions  of  x/a  and  y/a.  To  reveal  the  details  of  each  field  components 
clearly,  it  is  necessary  to  use  dilTercnt  scale  factors  for  different  plots.  To  show  the 
relative  amplitude  of  |e,(|  and  |cy|.  Ihe  ratio  of  |  peak/ 1  ®yl  peak 
I  ®yl  peak/1  ®x|  peak  indicated  in  the  plots.  For  the  lowest  order  polarization  mode 

(Figures  2a  and  2b),  the  peak  of  le„|  is  approximately  128  times  stronger  than  that  of 
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|ey| .  Ou  the  other  hand,  for  the  second  lowest  order  polarization  mode  of  the  same 
Gber  and  with  the  same  V  value  (F'igures  3a  and  3b),  |ey|  is  much  stronger  than  that 
of  |ej(|.  Comparison  of  these  plots  also  reveal  that  that  the  |ej(|  and  |  of  the 
lowest  polarization  mode  are  very  similar  to  that  of  |  Cy]  and  |  e,|  of  the  second  lowest 
order  polarization  mode.  Similar  results  have  also  been  obtained  for  the  same  fiber 
with  different  V  value  or  for  different  fibers.  In  summary,  we  note  for  each  polarization 
mode,  both  field  components  are  present,  one  field  component  is  much  stronger  than 
the  other.  In  addition,  the  dominant  field  component  ([exl  Figure  2a  and  |ey|  of 
Figure  3b)  has  a  peak  at  the  center  of  the  fiber,  just  like  those  of  the  LPjj  modes  for 
the  conventional  isotropic  fibers.  The  accompanying  and  much  weaker  field  com¬ 
ponent  (i.e.|ey|  of  Figure  2b  and  |e,j|  of  Figure  3a  respectively)  has  an  entirely 
different  look:  a  null  at  the  fiber  center,  surrounded  by  four  peaks,  one  in  each  qua¬ 
drant.  Our  calculations  also  show  that  fields  extend  much  more  into  the  cladding 
region  when  the  V  value  is  small,  i.e.,  the  fibers  are  operating  near  cutoff.  Of  course 
this  is  expected. 

For  the  fiber  discussed,  >  n,y  and  nj„  >  n^y.  and  the  dominant  electric  field 
of  the  lowest  polarization  mode  is  e„.  On  the  other  hand,  for  the  second  lowest  order 
polarization  mode,  the  dominant  electric  field  is  Cy.  In  general,  for  a  given  fiber,  the 
lowest  order  polarization  mode  is  the  mode  with  the  dominant  electric  fields  along  the 
direction  with  the  large  index  of  refraction.  Since  this  polarization  mode  has  the  larg¬ 
est  propagation  constant,  it  will  be  ideiP'fied  as  mode  s.  The  second  lowest  order 
polarization  mode,  labeled  as  mode  f,  has  its  dominant  electric  field  along  a  direction 
perpendicular  to  that  of  mode  s.  The  normalized  b  parameters,  as  defined  in  (Sla)  and 
(54b),  of  these  two  polarization  modes  are  labeled  as  bg  and  bf.  Curves  of  bg  vs  V  for 
various  values  of  d„y,.  and  d^yj  etc.  are  indi.stinguishable,  as  shown  in  f'igure  4.  llow'- 
ever,  curves  for  bf  do  vary  with  d,(y,  and  d„yj.  Of  particular  interest  is  the  bg  -  bf.  In 
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Figures  5,  6  and  7,  bg  -  b;  is  plotted  as  a  function  of  V  with  d^^^  and  as  parame¬ 
ters.  In  Figure  5,  djjy,  is  kept  to  a  constant  while  d^yj  varies.  We  note  that  near 
cutoff,  bg  -  bf  is  appreciable  if  djjyj  is  large.  This  is  understandable.  When  a  fiber  is 
operating  near  cutoff,  fields  spread  out  to  the  cladding  region,  and  therefore,  the  clad¬ 
ding  anisotropy  d„yj  should  and  do  have  an  strong  influence  on  the  fiber  characteristics. 
On  the  other  hand,  for  fibers  operating  far  from  cutoff,  fields  are  concentrated  in  the 
core  region  and  therefore  d^jy^  has  a  large  impact  on  the  bg  -  bj.  This  is  shown  in  Fig¬ 
ure  6  where  d,(yj  are  kept  constant  while  d^^y^  varies.  In  Figure  7,  d^^y,  =  d^yj  and  they 
vary  simultaneously.  For  these  cases  bg  -  b^  is  practically  independent  of  V. 


6.  Coupling  with  Biaxial  Fibers 

Although  the  excitation  of  low-birefringent  fibers  by  plane  waves  or  gaussian 
beams,  or  the  coupling  between  these  fibers  with  planar  or  channel  waveguides  have 
been  studied  extensively,  the  corresponding  problems  involving  high-birefringence  fibers 
have  not  been  treated  previously.  The  main  stumbling  block  is  the  lack  of  an  accurate 
description  of  fields  in  high  birefringence  fibers.  With  the  theory  developed  in  this 
report,  such  a  difficully  no  longer  exists  and  one  shall  be  able  to  proceed  to  study  these 
coupling  problems.  In  the  following,  we  shall  outline  the  basic  procedure  for  treating 
the  coupling  problems  and  present  some  of  the  results. 

Let  the  transverse  field  components  of  the  polarization  mode  u  be  and 
Then  the  transverse  electric  fields  in  the  fiber  can  bt  expressed  as 


F'or  a  given  incident  field  EJjp,  c„  is  given  by 
c.  =  f0"  X  K)  i,  ds 


(55) 


(56) 


where  a,  is  the  unit  vector  in  the  z  direction.  Thus  the  dependence  of  c^  as  a  function 
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of  the  state  of  polarization  (SOP)  of  the  incident  field  and  under  various  excitation 
conditions  can  be  determined.  We  assume  that  and  are  normalized  so  that  the 
power  coupled  to  the  polarization  mode  i/  is  proportional  to  |  c^| 

This  approach  has  been  employed  to  study  the  excitation  of  low  birefringent  as 
well  as  high  birefringent  fibers.  The  contrast  between  these  two  types  of  fibers  is  quite 
intere.sting.  Take  the  case  of  plane  wave  excitation  as  an  example.  As  the  input  SOP 
of  changes,  the  SOP  of  the  fields  in  the  low-birefringent  fibers  changes  with  the  input 
SOP,  but  the  total  power  coupled  into  the  fibers  remains  unchanged.  On  the  other 
hand,  for  high-birefringence  fibers,  both  the  SOP  and  the  total  power  vary  as  the  input 
SOP  varies. 
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Appendix  A 


In  this  Appendix,  the  expressions  of  foo(r)  etc.  are  listed 


I.  Core  region: 


Wr)  =  .io(^  u„  tI  +  T  IT  u„  f) 

■  II  -V 
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Appendix  B 

The  expressions  for  the  matrix  elements  of  Q  matrix  are  given  in  this  Appendix. 
Functions  foo(r)  etc.  are  listed  in  Appendix  A.  A  '  signifies  differentiation  with  respect 
to  r. 
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Figure  1.  Biaxial  fiber 


Figure  3.  Electric  Q 
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Figure  4.  Dispersion  of  biaxial  fibers 
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Figure?.  Birefringence  of  biaxial  fibers  with 


